Common origin of the two types of magnetic fluctuations in iron chalcogenides 
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We use inelastic neutron scattering to study the low energy spin excitations in moderately doped 
non-superconducting Fe1.01Teo.72Seo.2s- The spin excitations in this system contain components 
near (0.5,0,0) and (0.5,0.5,0) in a-b plane reciprocal lattice units using tetragonal unit cell notation 
(a=b=3.772 A and c=6.061 A). At low energies the scattering is centered around (0.5,0,0). With in- 
creasing energy, the spectral weight of low energy spin excitations centered around (0.5,0,0) abruptly 
shifts around 3 meV to the incommensurate spin excitations centered around (0.5,0.5,0). However 
both types of spin fluctuations exhibit the identical temperature dependence. These results indicate 
that the (0.5,0,0) type spin excitations and the incommensurate excitations around the (0.5,0.5,0) 
position have a common origin and both must be taken into account to understand the nature of 
magnetism and superconducting pairing in the iron chalcogenides. 

PACS numbers: 75.47.-m, 71.70.Ch 



I. INTRODUCTION 

The iron chalcogenides Fei+yTei-aSe^ have the a- 
PbO structure, which contains layers of Fe squares with 
the chalcogen atoms residing alternately above and below 
the centers of the squares, the same way in which FeAs 
layers are formed in the iron pnictides 1 -' 2 -. The supercon- 
ductivity arises when sufficient Se replaces Te in the an- 
tiferromagnetically ordered parent phase Fei+yTe 2 -. De- 
spite the similarities in crystal structure and Fermi sur- 
face topology 3 '-, the parent compounds of the iron pnic- 
tides and iron chalcogenides have very different mag- 
netic structures. The pnictides have a single stripe in- 
plane collinear (C-type, shown in Fig. 1(a)) antiferro- 
magnetic (AFM) structure^ characterized by the wave 
vector (0.5,0.5,0.5) in the notation of the tetragonal lat- 
tice as highlighted by the shaded area in Fig. 1(a), 
which coincides with the Fermi surface (FS) nesting wave 
vector— and the wave vectors of neutron spin resonance 
in superconducting samples*"*. The nonsuperconduct- 
ing Fei +y Te has a diagonal double-stripe bicollinear or- 
der (EMype, Fig. 1(b)) modulated along the (0.5,0,0.5) 
directio n 1 ^ 11 , whose in-plane component is 45° away 
from the FS nesting wave vector (0.5,0.5,0) where, cu- 
riously enough, the spin resonance of superconducting 
Fei+yTei-^Sea; is found 12 - 13 . The disparity between the 
static order and Fermi surface nesting as well as the large 
ordered momen t 11 ! 14 ' 15 in iron chalcogenides, have fu- 
eled the already heated debate about the nature of mag- 
netism in the iron based superconductors. The contro- 
versy over the magnetism in the iron chalcogenides has 
been mainly centered on whether it originates from itin- 
erant electrons^-^, localized momenta 1 ^— or both 2 ^— . 

As a good probe to magnetism, the spin dynamics 



in the iron chalcogenides has been extensively studied 
with inelastic neutron scattering (INS) measurements. 
In undoped Fei +y Te, the low energy excitations start ex- 
actly or closely to the (0.5,0,0) position depending on 
the amount of excess Fe 2 ^— . However, this dispersion 
about the magnetic zone center disappears at higher en- 
ergies, and new rings of excitations emerge above 60 meV 
around integer positions such as (1,0,0). The ring cen- 
tered at (1,0,0) disperses inward and eventually becomes 
a spot before disappearing above 275 meV— . In the 
27% Se-doped non-superconducting (NSC) sample 3 ^, the 
magnetic response around (1,0,0) starts from the lowest 
measured energy and forms incommensurate magnetic 
(ICM) quartets instead of a ring, but with increasing 
energy it evolves into a ring. The dispersion is still steep 
and persists to energies beyond 250 meV. The normal 
state magnetism in superconducting FeTei-^Sea^Sr— is 
very similar to that of the non-superconducting com- 
pounds except that the quartets never become a ring 
at high energies™ The observed a-b plane spin exci- 
tations in Fei+ y Tei- x Se x are summarized in Fig. 1(c) as 
schematics in reciprocal space. 

In most Se-doped compounds the (0.5,0,0) excitations 
are still present at low energies. Their intensity dimin- 
ishes with increasing Se-doping and becomes very weak in 
superconducting sample a 30 ' 34 ' 35 . It is therefore believed 
that the (0.5,0,0) type spin correlations have a deleteri- 
ous effect on the superconducting pairing. Since chemical 
inhomogeneity, impurity and phase separation exist in 
these materials^—, it is unclear whether these (0.5,0,0) 
spin fluctuations are intrinsic to the system or are a re- 
sult of an undesirable phase or domains segregated from 
the primary magnetism. 



To address the relationship between these two types 
of coexisting magnetic fluctuations, we have carried 
out INS measurements on moderately doped non- 
superconducting Fe1.01Teo.72Seo.28- Our results indicate 
that the (0.5,0,0) type excitations dominate the lowest 
energy spectral response, but then the strength of the 
scattering abruptly shifts to the ICM excitations cen- 
tered about (0.5,0.5,0), while the two types of spin fluc- 
tuations exhibit the identical temperature dependence. 
Taken together, we conclude that the two types of mag- 
netic excitations have a common origin. 

II. EXPERIMENTAL 

High quality single crystals of a-phase Fei+yTei-^Se^ 
were prepared with nominal composition of x — 0.3 us- 
ing the flux method and co-aligned with neutrons. The 
actual compositions of the samples were determined with 
prompt gamma activation analysis (PGAA) on beamlinc 
NG-7 at the NIST Center for Neutron Research (NCNR). 
Neutron scattering measurements were conducted on the 
cold neutron triple axis spectrometers (TAS) MACS and 
SPINS and the thermal TAS BT-7 at NCNR. Pyrolytic 
graphite (PG) was used as monochromator and analyzer 
for all the measurements and as filter for BT-7 measure- 
ments. A BeO filter was used on SPINS with fixed fi- 
nal energy Ef of 3.7 meV and horizontal collimations of 
open-80'-S-80'-open. Double focusing monochromator, 
Be filter, and all 20 channels of the detection system^ 
were employed on MACS with Ef fixed at 5 meV. Fixed 
Ef of 35 meV and Open-25'-S-25'-120' collimations were 
used for INS measurements on BT-7 4Q . Open-25'-S-25'- 
50' with Ef—1A.7 meV was used for elastic scattering on 
BT-7. The momentum transfer Q at (q x ,<ly,Qz) is de- 
fined as (H,K,L) = (q x a/2n, q y b/2w, q z b/2ir) reciprocal 
lattice units (r.l.u.) in the tetragonal unit cell {PA/nmm 
space group). The lattice parameters of the tetragonal 
unit cell are a = b = 3.772 A and c = 6.061 A at T = 1.5 
K. Elastic measurements were taken in both the (H, K, 0) 
and (H, 0, L) scattering planes, while the INS measure- 
ments were concentrated in the (H, K, 0) plane. The er- 
ror bars shown in the figures are statistical in nature and 
represent one standard deviation. 

III. RESULTS AND DISCUSSION 

We first determine the actual stoichiometry of our 
Fei+j,Tei_2;Se x crystals since these can differ significantly 
from the nominal compositions^. For this purpose, 
prompt gamma activation analysis (PGAA) was carried 
out on a small piece of single crystal. PGAA is a nonde- 
structive technique using neutron adsorption to simulta- 
neously determine the presence and accurate quantities 
of various elements in a compound 42,4 -. We find that 
the Fe:Te:Se molar ratio is 1.009:0.721:0.279. Relative 
expanded uncertainties for PGAA data are estimated at 
less than 5%. 

Neutron diffraction measurements were carried out on 
a small piece of single crystal from the same batch as the 
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FIG. 1: (Color online) (a) The schematic C-type collinear 
AFM order of the Fe moments in the a-b plane for iron pnic- 
tides such as LaFeAsO 4 ^. The shaded area represents the 
tetragonal unit cell which is used in the study, (b) The ab- 
plane projection of the E-type bicollinear AFM spin structure 
for iron telluride FeTe. (c) Schematics of the observed spin ex- 
citations for Se-doped FeTe in reciprocal space. The red star 
shows the position where FS nesting and the spin resonance 
has been observed for both pnictides and chalcogenides. The 
solid blue circles represents the low energy spin fluctuations at 
(0.5,0,0) and the filled green ellipses show the ICM excitations 
centered around the (0.5,0.5,0) position, (d) Temperature de- 
pendence of the integrated intensity of scans along (H,0,0.5) 
taken on BT-7. The inset shows the background-subtracted 
H-scans at some typical temperatures, (e) Constant-Q scans 
at (0.5,0.5,0) measured on SPINS at T=1.4 K and T=20 K. 
There is no evidence of a spin resonance or the development 
of a spin gap. 
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FIG. 2: (Color online) Contour plots of neutron scattering 
intensity in the (H,K,Q) plane at energy transfers of (a) 1.5 
meV (b) 3.5 meV, (c) 4.5 meV, (d) 6.0 meV (e) 7.0 meV and 
(f) 8.0 meV. The data were collected at T = 1.5 K on MACS. 
All the panels are plotted with the same color scale to show 
the intensity variations. The arrows in (a) and (b) show the 
directions of constant- E cuts plotted in Fig. 3. 
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crystals used for the inelastic neutron scattering. The 
crystal was aligned in the (H, 0, L) scattering plane. H- 
scans were performed at various temperatures. The static 
magnetic order survives in this compound as short range 
static order centered at the wave vector (0.47,0,0.5). This 
is consistent with the previous reports that Se-doped bulk 
non-superconducting samples have short range order cen- 
tered at incommensurate wave vector (0.5-<5,0,0.5) 34 i 44 ' 47 . 
The 5 value, 0.03 in this case, can be tuned by both the 
Se and Fe concentratio n 10 ! 44 , il-scans at some typical 
temperatures are shown in the inset of Fig. 1(d). The 
magnetic peaks are much broader than the instrumen- 
tal resolution with Lorentzian fits giving a width of 1.6 
(r.l.u.) corresponding to in-plane correlation length of 
3.7(8) A at 5 K. This is in agreement with that of sim- 
ilar doping FeTeo.ySeo.a^ 4 -. The integrated intensity of 
the ii-scans at different temperatures is plotted in Fig. 
1(d). Upon warming the peak intensity is gradually sup- 
pressed without an abrupt transition or change of peak 
position. The magnetic intensities cannot be detected 
above 60 K. The diffuse nature of the short-range or- 
der is reflected by the concave shape of the intensity- 
temperature curve. Constant- Q scans of inelastic neu- 
tron scattering at (0.5,0.5,0) at 1.4 K and 20 K are shown 
in Fig. 1(e). The overall spectrum shows little temper- 
ature dependence except some enhancement at 20 K at 
energies below 4 meV due to thermal population. Nei- 
ther a spin resonance nor the development of a spin gap 
at low temperature can be found, confirming the absence 
of bulk superconductivity. 

The above characterizations place our sample in the in- 
termediate doping part of the phase diagram where the 
long range AFM order is suppressed but bulk supercon- 
ductivity has not emerged. In this region, there are re- 
ports about weak charge carrier localization^ and spin 
glass ordering 4 ^"—. The magnetic and superconducting 
properties are sensitive not only to the Te/Se ratio, but 
also to the variation of Fe content^^^^ 4 -^ 1 -. The 
modest amount of excess Fe in our sample ensures that 
no significant complications arise due to the interstitial 
iron. 

Twenty grams of single crystals were co-aligned for the 
INS measurements which focused on the (H, K, 0) scat- 
tering plane because of the weak I/-dependenc o 27 i 30 ' 32 . 
For those measurements we drop the L coordinate for 
simplicity and present the data in terms of (H, K) only. 
In order to simultaneously monitor the two excitations 
we need to survey a wide region of momentum space in 
the low energy transfer range. MACS is ideal for this 
type of measurement because of the multiple detection 
systems—. An empty sample holder in the same sam- 
ple environment was also measured and used as back- 
ground subtraction. Fig. 2(a)-(f) show some typical con- 
tour plots at 1.5 K using data folded into the first quad- 
rant of the scattering plane. At low energy transfers the 
(0.5,0) type excitation is dominant with a broad peak. 
At the equivalent position of higher Brillouin zones, such 
as (1,0.5) and (1.5,0), peaks are weaker because of the 
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FIG. 3: (Color online) Contour plot of combined cuts for 
the T=1.5 K data in the (a) [0.5, K] direction and (b) the 
[H, 1 — H ] direction through (0.5,0.5) as a function of energy. 
Both figures are plotted on the same energy scale so that the 
correspondence between the two excitations can be seen. The 
open circle data in (b) are the result of a two-gaussian fit to 
the peaks shown in Fig. 5(a), which are measured using the 
high resolution configuration on SPINS. 
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FIG. 4: (Color online) (a) Constant-Q cut along E at (0.5,0) 
[black solid squares] and (1,0.5) [blue sold circles], (b) Cut 
at (0.5,0.5) along E. The shaded area is the energy range 
where the strong diffuse scattering is present. The dashed 
line shows roughly where the abrupt changes of intensities 
occur, (c) Integrated intensity of the cut through (0.5,0) along 
the H-direction with background subtracted, (d) The sum of 
fitted areas of the two peaks from the transverse scan through 
(0.5,0.5). Data obtained from BT-7 are used for E>8 meV. 



decreased magnetic form factor. As the energy increases 
the (0.5,0) type scattering quickly diminishes in intensity 
while the two ICM peaks around (0.5,0.5) start to appear 
and intensify. These two peaks are from separate sets of 
quartets of scattering about (1,0) and (0,1) respectively. 
We hereby call the excitation represented by these two 
peaks the ICM excitations to distinguish from the (0.5,0) 
excitations, and to avoid confusion about the Brillouin 
zone center. Also visible is the acoustic phonon mode 
stemming from (1,1) at low energies and from (2,0)/(0,2) 
at higher energies. 

At various energies, we performed cuts along [0.5, if] 
through the centers of both (0.5,0) and (0.5,0.5), which 
enables us to track the evolution of the two excita- 
tions simultaneously. Constant- E cuts in the transverse 
[H, 1 — H] direction through (0.5,0.5) were also conducted 
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which yield two symmetric ICM peaks that gives the 
dispersion of the ICM excitation. The combinations of 
these two cuts as a function of energy are shown in Figs. 
3(a) and 3(b). The [0.5, if] cuts at energies above 7 
meV are not included because the low-Q space could not 
be reached with the chosen fixed final energy. Fig. 3 
(a) shows that the disappearance of the broad peak at 
(0.5,0), as energy increases, is a three-stage process in- 
stead of a continuous one. The broad quasi-elastic scat- 
tering shows a sharp drop of intensity at about 1.5 meV, 
followed by another abrupt drop at about 3 meV be- 
fore the signal disappears above 6 meV. Fig. 3(a) also 
shows the gradually emerging excitation at (0.5,-0.5) and 
(0.5,0.5). The complete disappearance of the (0.5,0) exci- 
tation at around 6 meV is correlated with the opening of 
a spin gap for the ICM peaks around (0.5,0.5), as shown 
in Fig. 3 (b). In this partial gap the ICM peaks around 
(0.5,0.5) are suppressed in intensity, but their peak pro- 
files remain. Remarkably, the energies where the twin 
ICM peaks abruptly change spectral weight coincides 
with those where the opposite sudden change occurs for 
the (0.5,0) spectrum. This reciprocal interplay between 
the two spin excitations can be more easily seen in Fig. 4. 
The constant-Q scans at these two wave vectors show an 
opposite energy dependence. The presence of the (0.5,0) 
fluctuations is compensated by the gap opening in the 
ICM fluctuations near (0.5,0.5). The dashed lines in Fig. 
4(a) indicate the energies at which the two spectra show 
abrupt changes. The intensity of the constant-Q scan at 
(1,0.5) is also plotted on a log scale in Fig. 4(a) confirm- 
ing the energy dependence of the (0.5,0) correlation. 

In order to obtain the integrated intensity for the 
(0.5,0) spectrum, we performed cuts through (0.5,0) in 
the if-direction, instead of the if-direction to avoid the 
component of the excitation near (0.5,0.5). Because of 
the unreachable Q-space at higher energies, only half 
of the thereby obtained peak is integrated and plotted 
against energy in Fig. 4(c). The integrated intensity of 
the two ICM peaks around (0.5,0.5) is also plotted in Fig. 
4(d). The constant- E scans in the transverse direction 
through (0.5,0.5) were also performed with the thermal 
triple axis instrument BT-7. The raw data up to E = 
16 meV and lines of fits to two Gaussians are displayed 
in Fig. 5(b). The ICM excitation shows very little vari- 
ation in dispersion and in spectral weight between 8 and 
16 meV. The sum of fitted areas of the two Gaussians, to- 
gether with that of the constant- E cuts from the MACS 
data shown in Fig. 2, is plotted as a function of energy in 
Fig 4(d). Identical scans at E = 8 meV using the two in- 
struments were used to normalize the overall intensities. 
The spectrum of the ICM excitations about (0.5,0.5) re- 
mains constant above 7 meV. This is important because 
it means that the gap of the ICM excitations below 7 meV 
is not compensated by spectral gain at higher energies at 
the same wave vector as is the case in the superconduct- 
ing compounds. This again indicates that spectral weight 
is transfered between the two wavevectors and that the 
(0.5,0) spectrum is at cost of the ICM spectrum. Clearly, 
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FIG. 5: (Color online) Typical scans along the [H,l-H] direc- 
tion at various energy transfers carried out on (a) SPINS and 
(b) BT-7 at T — 1.5 K. The lines are Gaussian least squares 
fits convoluted with instrument resolution. The intensities 
and fits are shifted up sequently by an equal amount between 
adjacent scans for clarity. 



the spectral weight for these two types of excitations is 
inversely correlated, which rules out electronic phase sep- 
aration or magnetic inhomogeneity as the origin of the 
two types of magnetic correlations^. 

Another conspicuous feature in Fig. 3(b) is the disper- 
sion toward (0.5,0.5) below 3 meV before its steep out- 
ward dispersion. This hourglass shaped dispersion has 
been observed in both superconducting and over-doped 
non-superconducting system s 33 ' 50 . For this under-doped 
sample the inward dispersion is less pronounced and is 
more of a bell shape. To confirm this inward dispersion 
and the intense quasielastic scattering below 1.5 meV, we 
repeated the [H, 1 — H] scans on the cold neutron TAS 
instrument SPINS which, with Ef—3.7 meV and horizon- 
tally flat monochromator, offers better energy and Q res- 
olution. Fig. 5(a) shows these transverse scans through 
(0.5,0.5) at various fixed energy transfers. The low back- 
ground and distinct profile of the two peaks at E = 1 
meV in Fig. 5(a), in contrast to the broad quasielas- 
tic scattering extending to E — 1.5 meV in Fig. 3(b) 
and Fig. 4(b), is the result of improved resolution. The 
peaks are fit with two Gaussians convoluted with the 
instrumental resolution. The Q-positions obtained are 
over-plotted in Fig. 3(b) with open circles. The results 
are consistent with the MACS data and the bell shaped 
dispersion is clearly visible. The inward dispersion stops 
around 3 meV. Comparing our data with the studies of 
other compositions 3 - 3 -^, it seems the increasing Se dop- 
ing pushes the saddle point to higher energies. It should 
be also noted that in the energy range where the bell 
shape of the ICM excitation occurs, there is an abnormal 
change of the (0.5,0) spin spectrum as shown in Fig. 3(a). 
The change happens to both the line- width (inset of Fig. 
4(c)) and the peak intensity (Fig.4(a)). These anomalies 
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FIG. 6: (Color online) Temperature dependence of the cut 
(a) along [0.5,K] at E = 1 meV and (b) along [H,l-H] at 
E = 1 meV and (d) E = 10 meV. Note the change in the 
vertical color bar scales made to better show the change of 
intensities, (c) Temperature dependence of the full-width-at- 
half-maximum (FWHM) of the (0.5,0) spectrum at E = 1 
meV. 



are confirmed by a similar Q-E plot (not shown) for the 
cut through (1,0.5), which shows identical anomalies in 
intensity and linewidth. 

We now turn to the temperature dependence of the 
spin excitations. The (H, K, 0) planar maps have been 
obtained at different temperatures between 1.5 K and 
308 K for four typical energy transfers: 1 meV in the 
quasielastic region, 4.5 meV in the spin gap, 7.0 meV at 
the verge of the gap and 10 meV above the gap. Fig. 6 
shows the combined [0.5, K] cuts at E =1 meV, and the 
[H, 1 — H] cuts at E = 1 meV and 10 meV as a function 
of temperature. The integrated intensity of these cuts at 
all the above-mentioned energies are plotted in Fig. 7 
(a) and (b). Both types of magnetic correlations are so 
robust that they maintain their well-defined features up 
to the highest measured temperature for all the energy 
transfers. On warming the (0.5,0) spectrum at E = 1 
meV starts to gain intensity at about 60 K, where the 
static order disappears, reaches its maximum at about 
80 K and gradually decreases at higher temperatures. 
The ICM spectrum around (0.5,0.5), however, is clearly 
gapped below the transition temperature of the static 
AFM order, as shown in Fig. 6(b) and Fig. 7(b). As 
the system is heated from the short-range static ordered 
phase into the paramagnetic phase the static component 
also transfers to the background, resulting in the abrupt 
broadening of the linewidth of the (0.5,0) spectrum, as 
shown in Fig. 6(c). At E=7 meV the intensity of the 
ICM spectrum is less affected by the static order below 60 
K. At E = 10 meV (Fig. 6(d)), the ICM intensity remains 
unaffected by the static order and gradually increases 
monotonically all the way above 300 K. 

The neutron scattering intensity St{Q,w) at temper- 
ature T is related to the imaginary part of the dynamic 
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FIG. 7: (Color online) Temperature dependence of the inte- 
grated intensities of (a) the (0.5,K) scan at E = 1 meV and E 
= 4.5 meV, and (b) the [H, 1 - H] scan through (0.5,0.5) at 
E = 1 meV, 4.5 meV, 7 meV and 10 meV. (c) Temperature 
dependence of the dynamic susceptibility in the paramagnetic 
phase (60 K<T<310 K) obtained from the (0.5,0) type of ex- 
citations and (d) the (0.5,0.5) ICM excitations. The lines in 
(c) and (d) are fits to the power law x" = CT^ . 



susceptibility x(Q> w )t through St(Q,w) = [n(oj,T) + 
1]Xt(Q> w )i where n(u),T) is the Bose factor. In order to 
investigate the temperature dependence of the dynamic 
susceptibility in the paramagnetic phase, the St{uj) in- 
tensities for T> 80 K in Fig. 7(a) and (b) are converted 
to the values of dynamic susceptibility xi^Yr an< ^ plotted 
in Fig. 7(c) and (d). The scattering at energy transfers of 
7 meV and 10 meV increases with increasing temperature 
as expected due to the thermal factor, but the susceptibil- 
ity decreases with T indicating that the intrinsic strength 
of the magnetic scattering decreases as the scattering 
evolves to higher temperature. The T-dependencies of 
the dynamic susceptibility in the paramagnetic phase 
(T>60) were fit to the power law X " = CT P . The (3 
values obtained for wavevectors (0.5,0) and (0.5,0.5) at 
E=l meV, denoted by black squares in Fig. 7(a) and 
7(b) (solid and empty), are -1.28±0.02 and -1.26±0.03 re- 
spectively. Similarly the /3 values at i?=4.5 meV (red cir- 
cles) for these two vectors are -0.68±0.01 and -0.65±0.01. 
Note that the dynamic susceptibility of the two excita- 
tions have the identical temperature dependence. This 
is further evidence that both types of excitations have a 
common origin. At an energy of 7 meV (Fig. 7(b)) where 
the (0.5,0) spectrum is almost completely depleted, the 
x" for (0.5,0.5) continues to follow the same relationship 
with /3=-0.66±0.01. At an energy transfer of 10 meV, /3 
becomes -0.54±0.03. 
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For the superconductors, the high energy magnetic ex- 
citations remain unchanged when cooled from the normal 
state to below T c 32 ' 50 . It is the low energy part of the 
magnetic excitation spectrum that responds to the for- 
mation of superconducting pairs, as would be expected. 
These changes include the opening of a spin gap and the 
development of a spin resonance. In our NSC system, the 
gap still develops for the ICM spectrum centered around 
the same (0.5,0.5) position, but without the development 
of superconductivity or a spin resonance. Instead, we 
have the (0.5,0) type excitations, which apparently cor- 
respond to that part of magnetism that is needed for the 
spin resonance in the superconducting state, as only the 
(0.5,0) type excitations are suppressed when the super- 
conductivity and the associated spin resonance develop. 
Now that we know that the two type of excitations have a 
common origin, suppressing the (0.5,0) type correlations 
through doping Se cannot be simply understood as elim- 
inating a coexisting phase. Rather one has to treat the 
two excitations as one problem when trying to reveal the 
driving force for magnetism and superconductivity. In 
that sense, the itinerant electrons alone may not be able 
to provide a complete answer. The inter-band Fermi sur- 
face nesting describes the main features of the ICM exci- 
tations such as the incommensurate excitations 32 and the 
hour-glass dispersion near (0.5,0.5)^2, but the Fermi sur- 
face near the X point has not been founc&S 3 - yet to sup- 
port the nesting scenario for the (0.5,0) excitatio n 16 i 17 . 

In a local moment picture, the magnetic ground state 
is governed by superexchange interactions. The contest 



between the collinear and bi-collinear order is controlled 
by the competition between Ji, J 2 and J3, which have 
different chalcogen height dependencies^—. As Se re- 
places Te, the chalcogen height is reduced^ which results 
in increased Ji, J2, decreased J3, and consequently a less 
favored bicollinear order at (0.5,0). The spin wave spec- 
trum centered at (0.5,0) that extends up to 60 meV— in 
undoped FeTe is suppressed to lower energies, and is com- 
pletely taken over by the ICM magnetism as the Se con- 
tent increases well into the superconducting region. This 
scenario is supported by the reciprocal interplay between 
the two types of excitations and the temperature depen- 
dence of the \ow-E dynamic susceptibility presented in 
this study. However, it is not possible to understand the 
low energy features such as the hour-glass dispersio n 26 ! 33 
and the abnormal change of intensity and linewidth for 
the (0.5,0) spectrum with just a local spin picture. Our 
results call for a more unified mechanism that reconciles 
these features and embraces the two types of excitations 
as having a common origin. 
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